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Summary 


This  research  program  centered  on  the  investigation  of  low-accuracy 
analog  opto-electronic  methods  for  processing  multi-sensor-array  data  (e.g., 
phased  array  radar,  sonar).  Two  methods  were  considered  for  signal 
orthogonalization  in  low-accuracy  preprocessors,  one  based  on  a  resonant 
piezo-electro-optic  modulator,  the  other  based  on  photo-refractive  and 
similar  spatial  light  modulator  devices.  Neither  method  appears  to  offer 
significant  benefits  over  all-electronic  preprocessor  methods  because  of  the 
lack  of  large-scale  parallelism.  An  analog,  continuous-time  opto-electronic 
processor  for  adaptive  phased  array  signal  processing  was  also  considered 
that  does  exhibit  considerable  parallelism  and  shows  promise  for  use  as  a 
low-accuracy  co-processor  in  a  hybrid  analog/digital  configuration. 


S0WAR7  or  RKSBARCH 

1.  Objectives  of  Research 

The  principal  objective  of  this  research  program  was  to  determine  whether 
low-accuracy  analog  opto-electronic  signal  processing  methods  could  be  used 
effectively  to  pre-process  multi-sensor  array  signal  data  to  improve  the 
overall  performance  of  an  analog  optical/digital  electronic  algebraic  signal 
processing  system. 

2.  Background 

The  Air  Force  has  many  signal  processing  tasks,  particularly  those 
associated  with  multi-sensor-array  data,  that  lie  beyond  the  current 
capabilities  of  electronic  signal  processing  systems.  The  development  of 
highly  parallel  processing  architectures  and  algorithms,  coupled  with 
improvements  in  VHSIC  technologies,  will  increase  significantly  the 
capabilities  of  electronics  for  such  tasks.  Optics,  with  its  high  degree  of 
parallelism,  enormous  bandwidth,  and  three-dimensional  interconnect 
capability,  also  offers  potential  solutions. 

Emphasis  on  optical  signal  processing  research  has  been  placed  in  recent 
years  on  the  study  of  algebraically-oriented  processors  suited  for  specific 
operations  like  direction  finding,  beam  forming,  and  adaptive  noise 
cancelling,  and  for  more  general  operations  like  least-squares  solution  of 
sets  of  linear  equations  and  singular-value  decomposition.  Often,  high 
accuracy  (e.g.,  32  or  more  bits  on  a  digital  processor)  is  required  for 
satisfactory  performance.  Unfortunately,  optical  processors  are  inherently 
analog  in  nature,  and  digital  accuracy  is  obtained  only  at  a  significant  price 


in  terms  of  system  throughput,  complexity,  and  cost.  For  example,  the  SAOBiC 
processor,  at  one  time  under  development  by  GuilTech  Research  Company  [1], 
though  offering  in  excess  of  500x10®  32-bit  floating-point  operations  per 
second,  was  still  slow  compared  to  all-electronic  processors  and  far  behind 
the’  lO^-IO11  operations  per  second  offered  by  analog-accuracy  optical 
processors. 

One  well-established  method  for  increasing  the  speed  of  processing 
systems  for  data  from  multi-sensor  arrays  is  through  suitable  preconditioning 
of  the  input  signal  data.  Preprocessing  is  regularly  performed  on  signal  data 
in  adaptive  phased  array  radar  signal  processors.  Significant  increases  in 
overall  processor  speed  generally  result  and,  in  some  cases  at  least, 
reduction  in  system  cost  as  well  [2].  The  key  to  increasing  convergence  rate 
in  these  cases  is  in  effect  a  transformation  applied  to  the  input  signal 
correlation  matrix  that  reduces  the  spread  in  the  eigenvalues  of  the  matrix. 
This  is  equivalent  to  improving  the  condition  number  of  the  matrix,  with 
greater  stability  resulting  for  subsequently-applied  algorithms.  A 
particularly  important  aspect  of  such  preprocessing  is  that  it  need  not  be 
performed  with  high  accuracy:  low-accuracy  analog  preprocessors,  effecting, 
for  example,  an  approximate  orthogonalization  on  the  input  signal  vectors,  are 
regularly  used  in  cascade  for  adaptive  phased  array  radar  signal  processing 
[3] . 

Two  major  areas  of  investigation  received  emphasis  in  the  research 
program.  One  centered  on  the  basic  operations  required  of  opto-electronic 
components  for  signal  orthogonalization  and  on  ways  for  performing  those 
operations.  The  other  centered  on  a  particular  processor  for  performing 
adaptive  phased  array  radar  signal  processing.  Both  areas  are  discussed  in 
detail  in  the  following  sections. 

r 


3.  Signal  Orthogonal ixat ion 

One  of  the  most  attractive  preprocessors  for  adaptive  phased  array  signal 
processing  is  the  Gram-Schmidt  cascade  preprocessor  [4],  which  effectively 
decorrelates,  or  orthogonalizes,  signals  from  the  antenna  array  elements. 
Figure  1  shows  a  system  using  the  Howells-Applebaum  form  of  adaptive-loop 
processor  for  signal  orthogonalization.  A  single  adaptive  loop  is  shown  in 
Fig.  2.  The  purpose  of  the  loop  is  to  compute  an  output  signal  y^(t)  that 
contains  only  information  that  is  not  contained  in  signal  y^(t),  i.e., 

innovative  information.  This  is  accomplished  by  calculating  Yj(t>  in  accord 
with  the  equation 


Yj  (t)  -  X^t)  -  Y^ttjy.tt) 


where  (t)  is  a  correlation  coefficient  given  by 


<yi*(t)x.(t)> 


Yij(t)  "  <yi*(t)yj(t)> 


the  angle  brackets  denoting  a  time  average.  It  is  easily  shown  that  in  the 
steady  state,  yi(t)  and  y^(t)  are  uncorrelated,  i.e.,  <y^* (t)y^ (t) >  ■  0.  For 
the  case  of  interest,  y^(t),  x^(t),  etc.,  represent  the  complex  envelopes  of 
associated  narrowband  signals. 

It  is  evident  from  Eq.  (1)  that  an  opto-electronic  preprocessor  of  the 
Gram-Schmidt  type  must  be  able  to  perform  three  elementary  operations:  1) 
signal  multiplication;  2)  correlation  coefficient  calculation,  which  involves 
a  time  average;  and  3)  signal  subtraction. 


Two  basic  approaches  to  performing  these  three  operations  were  considered 
in  this  program.  One  involves  the  use  of  a  resonant  piezo-electro-optic 


Modulator  device.  The  other  relies  on  time-averaging  characteristics  of 
quasi-real-time,  light-addressable  spatial  light  modulators.  These  two 
approaches  are  discussed  in  the  following  sections. 

3.1.  Besonant  Pieso-electro-optic  Modulator 

A  resonant  piezo-electro-optic  modulator  [5,6,7]  can  perform  both  the 
multiplication  of  y^(t)  by  Y^ft)  in  Eq.  (1)  and  the  time  average  needed  to 
calcualte  Y^ft). 

The  device  is  similar  in  many  respects  to  an  ordinary  linear  electro¬ 
optic  modulator,  but  it  is  operated  in  an  a.c.  mode  at  the  piezo-electric 
resonance  vibration  frequency  of  the  crystal.  At  this  frequency,  shear  and 
thickness  mode  vibrations  of  the  crystal  result  in  large  sinusoidal 
fluctuations  of  the  birefringence  of  the  crystal  with  the  application  of 
relatively  small  a.c.  voltages.  If  placed  between  crossed  polarizers  and 
driven  at  the  resonance  frequency  ,  the  modulator  produces  fluctuations  in 
light  intensity  transmittance  t  in  accord  with  the  equation 


2  1 

sin  [  —  *  cos  u  t  ] 

2  o 


where  ♦,  the  peak  birefringent  retardation,  is  directly  proportional  to  the 
amplitude  of  the  applied  sinusoidal  voltage.  Depending  on  the  crystaline 
material  used,  resonant  frequencies  in  the  range  0.5  to  500  MHz  appear  to  be 
easily  obtained  [5,6].  Because  of  the  mechanical  resonance,  voltages  required 
to  achieve  significant  changes  in  transmittance  t  are  typically  several 
orders  of  magnitude  lower  than  those  required  for  the  d.c.  linear  electro¬ 
optic  (Pockels)  effect.  Specifically,  a.c.  voltages  in  the  5-20  volt  range 
are  adequate  in  reported  cases  for  high  efficiency  modulation. 


For  optical  processing  applications  device  transmittance  t  should  vary 


linearly  with  the  applied  signal.  If  the  crystal  used  has  some  natural 
birefringence  or  is  illuminated  by  light  of  the  appropriate  elliptical 
polarization,  the  operating  point  of  the  device  can  be  shifted  such  that 
t  has  the  form 


.  2  .  it  4> 
T  ■  sin  [  —  +  — 


cos  to  t  J 
o 


(4) 


Theoretical  considerations  show  that  for  an  applied  driving  a.c.  voltage  V  at 
an  arbitrary  frequency  to,  device  response  can  be  wr itten  as 


r  ■  sin2  [  ~  8  (  ~  )  cos  tot  ]  (5) 

o 

where  VQ  is  a  constant  that  appears  typically  to  be  in  the  1-20  volt  range  and 
8  is  a  resonance  parameter  given  (in  the  vicinity  of  resonance  ~  far  off 
resonance  t  is  essentially  constant)  by  the  standard  resonance  function 

1 

to-to  ,  - 

8  *  8  (to)  -  (  1  +  4QZ(  -j—  )  ]  (6) 

o 

Device  Q,  unless  deliberately  reduced,  is  typically  of  the  order  of  10^,  and 
the  reasonance  is  thus  extremely  sharp. 

So  long  as  the  applied  voltage  V  is  kept  sufficiently  small,  t  can  be 
written  in  the  form 

t  *  —  +  8  (to)  (  —  )  cos  tot  (7) 

o 


There  is  a  phase  shift  introduced  for  to  *  toQ 
results  in  the  applications  of  interest. 


r 


but  this  does  not  affect  the 


The  high-Q  resonance  just  described  can  be  exploited  in  estimating  the 
correlation  coefficient  y  of  Eq.  (2)  through  the  approximation 

Yij<fc)  “  T  f  y i* <5>y j  < 5)  exp  [  ld5  (8) 

Figure  3(a)  illustrates  estimator  operation  schematically  for  a  single  pair  of 
antenna  element  signals;  Fig.  3(b)  shows  the  same  operation  in  complex  signal 
form.  In  the  latter,  the  resonant  filter  is  represented  by  the  transfer 
function  associated  with  its  baseband  equivalent,  which  is  a  simple  RC 
filter.  The  output  of  the  filter  is  slowly  time-varying,  with  approximate 
bandwidth  1/T. 

It  should  be  emphasized  that  the  high-Q  filter  of  Fig.  3(a)  is  the 
electro-optic  device  itself  and  that  its  output  is  not  an  electrical  signal 
but  rather  the  modulated  light  transmittance  r  of  Eq.  (7).  Figure  4  models 
the  essential  input-output  characteristics  in  terms  of  a  resonant  tuned 
circuit  and  the  electo-optic  transfer  characteristics.  The  input  signal  is 
the  product  of  the  two  IF  signals  of  Fig.  3,  with  the  filter  tuned  to  the 
difference  frequency  uiQ. 

Figure  5  shows  a  possible  conf iguration  of  the  adaptive  loop  processor 
using  the  resonant  piezo-electro-optic  modulator,  with  all  signals  represented 
for  convenience  in  complex  form.  The  function  {b  +  y ,exp[  jcu^t] }  represents 
the  light  intensity  illuminating  the  resonant  piezo-electro-optic  modulator, 
bias  b  being  added  to  avoid  clipping  on  negative  parts  of  the  signal. 

The  overall  Gram-Schmidt  signal  orthogonalizer  would  consist  of  a  number 
of  such  loops  operating  in  a  configuration  suggested  by  Fig.  1.  System 
dynamic  range  is  limited  primarily  by  nonlinear  terms  not  contained  in  Eq.  (7) 
and  by  noise  associated  with  the  bias. 


3.2.  Photo-Addressable  Spatial  Light  Modulator 

As  noted  earlier,  the  resonant  piezo-electro-optic  modulator  serves  both 


to  calculate  the  correlation  coefficient  through  its  high-Q  time-averaging 
characteristics  and  to  multiply  (t)  by  that  coefficient  through  its  a.c. 
light  intensity  transmittance.  An  alternative  approach  is  to  use 
photorefractive  or  other  light-addressable  spatial  light  modulator  (SLM) 
mater ials. 

The  idea  is  based  on  the  observation  that  Young's  interference  fringes 
can  be  used  a  direct  measure  of  the  correlation  of  two  optical  waves.  In 
particular,  if  two  tilted  plane  waves  with  complex  amplitudes  y^(t)  and  y^(t) 
are  allowed  to  interfere,  the  resultant  fringe  pattern  has  amplitude 
proportional  to  the  correlation  coefficient  Y^  of  Eq.  (2).  If  this  fringe 
pattern  is  used  to  expose  a  light-addressable  SLM,  a  grating  structure  results 
that  can  diffract  an  incident  light  wave.  Assuming  linearity,  the  amplitude 
of  the  diffracted  wave  is  proportional  to  the  grating  amplitude,  which,  as 
noted,  is  proportional  to  Y  .  Thus  a  wave  with  amplitude  proportional  to 
Y^i  j  { t )  can  be  produced. 

One  possibility  is  to  use  a  pair  of  acousto-optic  modulators  driven  by 
narrowband  signals  having  complex  envelopes  y^(t)  and  y^ (t) •  These  produce  a 
pair  of  waves,  incident  on  a  light-addressable  SLM,  with  complex  amplitudes 

Ui(x,y,t)  -  Yi(t)  exp[j2irfQxJ  expfj^tj 


(9a) 


(9b) 


t^(x,y,t)  -  Yj(t)  exp[j2iffox]  explju>ot] 

The -interference  of  these  two  waves  produces  an  exposing  intensity  at  the  SLM 
of  the  form 


I(x,y)  a  bias  +  Iy^I  coslu^x  +  argtY^)] 


where  y^  is  the  correlation  coefficient  specific  to  a  time  integration  over 
the  exposure  period. 

Assuming  a  linear  response  of  the  SLM  to  the  exposing  wave  intensity,  a 
complex  wave  amplitude  transmittance  of  the  form 


t (x,y)  ■  t.  +  o(  y  expl j2irf  x]  +  Y*  exp[-j2*f  x]  ]  (11) 

DO  O 

results.  If  the  grating  is  illuminated  with  light  of  complex  amplitude 
y^ (t)exp[-j2'«rf0x] ,  a  transmitted  wave  proportional  to  y^(t)y^j  results, 
which,  with  a  proper  phase  relationship,  can  be  subtracted  from  a  wave  of 
amplitude  Xj(t). 

In  Eq.  (11)  y  is  written  as  a  constant.  In  general,  however,  y  is  a 
function  of  time,  varying  with  a  bandwidth  that  may  exceed  tens  or  even 
hundreds  of  kilohertz.  In  order  for  the  grating  transmittance  to  vary  with 
time  in  the  proper  way,  the  light-addressable  SLM  grating  must  "track"  the 
local  time  average  of  the  exposure  pattern.  Specif ically,  it  should  have  a 
time  constant,  both  for  fringe  buildup  and  decay,  in  the  range  from  tens  of 
milliseconds  to  tens  of  microseconds,  depending  on  the  application. 
Photorefr active  materials  such  as  BSO  exhibit  such  characteristics  and  would 


be  a  possible  candidate.  For  photorefr active  materials,  the  rise  time  £or 
writing  differs  from  the  relaxation  time  for  erasure  and  is,  in  fact, 
dependent  on  the  intensity  of  the  writing  light  distribution,  suggesting  the 
need  for  automatic  gain  control. 

Although  the  correlation  coefficient  and  the  product  yy^  can  in  theory 
be  calculated,  this  scheme  presents  too  many  difficulties  to  be  particularly 
practical  for  real-time  signal  processing  tasks.  For  one  thing,  the  wave 
field  being  diffracted  by  the  grating  must  not  seriously  perturb  the  grating 
structure  itself.  This  means  that  the  diffracted  wave  field  should  be 
incoherent  with  the  two  interfering  waves,  preferably  of  a  wavelength  at  which 
the  SLM  is  not  sensitive.  Probably  more  serious  a  limitation  is  that,  unlike 
for  the  resonant  piezo-electro-optic  modulator,  complex  wave  amplitudes  bear 
the  signal  information.  The  adaptive  loop  processor  must  therefore  have  the 
form  of  an  interferometer,  capable  of  modulating,  adding,  and  subtracting  wave 
amplitudes.  When  this  limitation  is  added  to  the  architectural  difficulties 
discussed,  the  SLM  scheme  appears  to  be  impractical  for  the  processing 
operation  of  interest  here. 

3.3.  Parallelism  Coos ider at ions 

In  the  final  analysis,  although  both  the  resonant  piezo-electro-optic 
modulator  scheme  and  the  light-addressable  SIM  scheme  could  perform  the 
desired  functions,  neither  appears  to  be  competitive  with  all-electronic 
schemes.  This  conclusion  is  reached  largely  through  architectural 
considerations.  Simply  put,  that  fraction  of  an  opto-electronic  Gram-Schmidt 
signal  orthogonal izer  that  would  be  optical  is  relatively  minor  and,  of  much 
greater  importance,  it  does  not  exploit  optical  parallelism  in  any  way.  A 
conclusion  reached  by  many  optical  signal  processing  researchers  is  that  an 
optical  processor  will  not  compete  with  its  electronic  counterpart  unless  the 


optical  processor  operates  with  a  high  degree  of  parallelism,  preferably 
through  global  interconnections.  The  Gram-Schmidt  processor  requires  none  of 
this  and  there  is  no  reason,  in  fact,  that  the  integration-multiplication 
capabilities  of  the  opto-electronic  devices  discussed  should  not  be  replaced 
by  electornic  multipliers  and  integrators. 

4.  Adaptive  Beamforming  Processing 

One  way  in  which  the  resonant  piezo-electro-optic  modulator  device  can  be 
used  with  much  more  optical  parallelism  is  in  the  direct  processing  of  antenna 
element  signals  for  adaptive  beamforming.  There  is  no  longer  the  sought-for 
relationship  between  a  low-accuracy  analog  preprocessor  and  a  high-accuracy 
digital  post-processor.  However,  it  may  be  possible  to  use  the  low-accuracy 
optical  system  in  a  co-processor  configuration  along  the  lines  recently 
proposed  by  Caulfield  [10]. 

The  system  for  performing  the  radar  signal  processing  calculations  is  an 
electro-optical  matrix-vector  multiplier  in  an  analog  feedback  system. 
Related  systems  have  been  investigated  before  [8] ,  but  always  they  have 
implemented  the  necessary  complex  arithmetic  using  relatively  awkward 
negative-real  representations  for  complex  variables.  Furthermore,  they  have 
operated  on  an  iterative  basis  rather  than  the  continuous  basis  facilitated  by 
the  resonant  piezo-electro-optic  modulator. 

The  principal  goal  of  the  processor  is  the  calculation  of  the  optimum 
complex  array  weighting  vector  w  that  satisfies  the  matrix-vector  equation 


s*  -  Mw 


(12) 


where  M  is  the  covariance  matrix  associated  with  the  (complex)  signal  from  the 
antenna  elements  and  s  is  the  array  steering  vector.  Iterative  processors 
calculate  w  via  the  algorithm 

*k.i  •  **  ♦  « -  «”> 

the  complex  arithmetic  being  implemented  using  nonnegative-real  (e.g.,  three- 
component  or  biased  real-imaginary)  representations  for  the  matrix  and 
vectors.  Figure  6  illustrates  a  system  for  calculating  w  using  dynamic 
relaxation  methods.  The  integrators  in  the  system  can  be  replaced  by  simple 
RC  filters  if  the  matrix  M  (which  represents  the  signal  environment  in  the 
adaptive  phased  array  radar  case)  is  time-varying.  It  can  be  shown  that  the 
system  is  stable  if  the  matrix  M  is  positive-definite,  in  which  case  the 
eigenvalues  of  M  have  real  parts  that  are  positive.  The  covariance  matrix  M 
of  Eq.  (12)  is  indeed  positive-definite,  and  non-oscillatory  convergence  is 
assured. 

Operation  of  the  system  is  described  in  terms  of  the  inversion  of  the 
matrix  equation 

y  ■  Mx  (14) 

where  x  is  to  be  solved  for.  In  the  radar  case,  all  vector  and  matrix 
components  are  narrowband  signals  and  can  be  written  in  the  form 


'rj'*’:  -r,  mv  irj  r;  ^  ^  '.■  mi  *w 


lx  (t)l  cos  (w  t  +  0  (t)  ]  (15) 

&  ox 

|yi(t)|  coa(uot  +  9y(t)J  (16) 

lm(t)|  co«t« t  +  e..(t)J  (17) 

o  13 

The  key  to  system  operation  is  an  SLM  made  up  of  resonant  piezo-electro-optic 
modulator  elements  exhibiting  the  time-varying  light  intensity  transmittance 

Tj.  ■  b  +  Irn^l  cos(uot  ♦  9^1  (18) 

where  b  is  a  bias.  Overall  system  operation  is  indicated  by  the  block  diagram 
of  Pig.  7.  Certain  mixer-filter  combinations  have  been  made  unnecessary  by 
the  use  of  complex  representations  for  the  signals.  The  optical  matrix-vector 
multiplier  is  in  most  respects  like  the  system  described  by  Goodman  et  al.  [9] 
except  that  the  matrix  elements  m^  are  conveyed  on  a  temporal  frequency 
carrier. 

The  SLM  representing  M  is  made  up  of  a  2-D  matrix  of  individual  resonant 
piezo-electro-optic  modulators  of  the  type  discussed  earlier.  The  input 
signal  to  the  ijth  element  is  the  product  of  the  ith  and  jth  IP  signals: 

V.^(t)  -  !xi (t)  I  cos{(w1  +  u>o)t  +  ©i ( t)  1  (19) 

|x^(t)  I  costo^t  +  9^(t)]  . 

with  the  filter  being  tuned  to  the  difference  frequency  .  Processing 


requires  that  all  signal  pair  products  be  calculated  and  applied  to  separate 


electro-optic  elements  of  the  matrix  SLM. 

Assuming  the  small  signal  approximation  leading  to  Eq.  (7)  is  satisfied, 
the  transmittance  of  the  ijfc**  cell  of  the  resonant  SIX  is  given  by  Eq.  (18). 
The  transmittance  function  evolves  with  a  bandwidth  approximated  by  1/Tr 
whereas  the  relaxation  algorithm  converges  at  a  rate  governed  by  the  smallest 
eigenvalue  of  M. 

S.  General  Conclusions 

Low  accuracy  opto-electronic  systems  can  be  used  for  preprocessing  of 
multi-sensor-array  data  to  reduce  the  load  on  high-accuracy  post-processors. 
However,  based  on  architectural  considerations  of  the  Gram-Schmidt 
preprocessor ,  it  would  appear  that  there  is  insufficient  parallelism  in  the 
optical  rystems  to  warrant  replacing  all-electronic  preprocessors.  The 
resonant  piezo-electro-optic  modulator  approach  considered  shows  promise  for 
other  signal  processing  applications  where  complex  arithmetic  must  be 
performed.  In  particular,  it  shows  promise  for  adaptive  beamforming  in  cases 
where  high  accuracy  is  not  required  or  where  analog  and  digital  coprocessors 
work  side  by  side. 
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